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We present results from a lattice study of SU (2) color, — 1 supersymmetric Yang-Mills theory 
using domain wall fermions. Supersymmetry in this particular lattice formulation is expected 
to emerge in the continuum and chiral limits without any fine-tuning of operators. Preliminary 
results for the static quark potential, residual mass, chiral condensate and spectrum-a potential 
indicator of supersymmetry restoration-are presented and discussed. 
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1. Introduction 

In recent years, much effort has been devoted toward formulating supersymmetric (SUSY) 
gauge theories on the lattice. This has been partially motivated by the theoretical as well as techni- 
cal challenges associated with the problem, as well as the obvious potential role of SUSY in beyond 
the standard model physics. Of crucial importance pertaining to the latter point is an understanding 
of dynamical symmetry breaking-something which may in principle be achieved with numerical 
simulations provided an appropriate lattice discretization of the theory is found. 

Since naive discretizations typically break SUSY explicitly, lattice simulations require fine- 
tuning in order to cancel off the undesirable SUSY breaking operators which may arise through 
radiative corrections. However, it has been realized for some time that one of the simplest SUSY 
theories, jV = 1 supersymmetric Yang-MiUs (SYM), may be simulated with conventional lattice 
discretizations and yet require only a minimal degree of fine-tuning. Although this theory does not 
possess SUSY breaking (as implied by a nonvanishing Witten index [^), it is believed to exhibit a 
variety of other interesting features such as discrete chiral symmetry breaking. The field content of 
^ = 1 SYM consists of a vector field and a single adjoint representation Majorana fermion, and in 
conventional lattice discretizations of this theory, the only relevant SUSY violating operator which 
may arise radiatively is a gluino mass term. As a result, in the chiral and continuum limits SUSY 
is restored accidentally at infinite space-time volume. 

In the past, a variety of numerical studies have employed Wilson fermions in order to simulate 
c/K = 1 SYM (for a review, see [^), however, these were subject to both difficulties of fine-tuning 
and the sign problem. In contrast, it was observed in [||] that domain-wall fermions (DWFs) are an 
ideal fermion discretization for simulating = 1 SYM because of their good chiral properties [Q, 
|5|, |6|] and positivity of the fermion Pfaffian obtained from "integrating out" the fermion degrees of 
freedom. In this formulation, the chiral limit may be achieved without any fine-tuning of operators. 

The first and until recently the only study to use DWFs to investigate .yV = 1 SYM focused 
on the chiral limit of the gluino condensate [^. In our work, we expand on this early study in 
several important respects: we 1) establish the lattice scale by measuring the static quark potential 
and provide evidence for confinement which is consistent with expectations, 2) determine the size 
of the residual mass in order to ascertain the proximity to the SUSY point, 3) extrapolate the 
chiral condensate to the chiral Umit using a recent, theoretically motivated fit formula for its 
dependence, and 4) study the spectrum of the theory. With exception to the third point, questions 
such as these could not be easily addressed in [0] because space-time volumes where too small. 

Finally, we note that a similar but independent study of =yK = 1 SYM was presented by J. Giedt 
at this conference [^. 

2. Simulation and measurement details 

Numerical simulations of ^ = 1 SYM were performed using an appropriately modified ver- 
sion of the Columbia Physics System (CPS). We use a Wilson gauge action with Majorana DWFs 
in the adjoint representation of the gauge group SU{2). The specific details of this lattice action 
may be found in [^. Rational Hybrid Monte Carlo simulations were performed on a 16^ x 32 x 
lattice with Lg = 16,20,24 and 28, gluino masses mf = 0.01,0.02 and 0.04, a domain-wall height 
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Figure 1: Monte Carlo time history for the gluino 
condensate for j3 = 2.3 and mf — 0.02. 



Figure 2: Heavy quark potential as a function of dis- 
tance for mf = 0.02 and — 16. Sommer scale error 
bars are indicated by dashed lines. 



of M = 1.9, and coupling jS = 2.3. Several additional simulations where performed at the weaker 
couplings jS = 2.3533 and j8 = 2.4 as well. For each ensemble, a total of 2500 to 3000 trajecto- 
ries were generated starting from an "ordered" configuration and equilibrium was achieved within 
the first 500 trajectories. Measurements were made using uncorrelated configurations which were 
generated thereafter. A plot of the gluino condensate time history is shown in Fig. (jl]). 



3. Results 

3.1 Static quark potential 

The static quark potential was extracted from Wilson loop measurements for a range of cou- 
plings. Wilson loops were measured in the fundamental representation of the gauge group and on 
Coulomb gauge fixed gauge configurations. For fixed distances r = |x|, the Wilson loops were fit 
as a function of time to the formula: 

(W(x,?)>=C(x>-^(''>, v{x)=Vo-^^+a\x\ (3.1) 

within an interval were excited state contamination appeared to be negligible. The extracted values 
of y (x) were then fit to the Cornell potential. The constant term (Vq), Liischer term (a), string 
tension (a), and Sommer scale (ro) defined by: 



1.65 (3.2) 



were determined by double jackknife fits to the data. Table (|I|) and Fig. (^ summarizes the fit 
results. A decrease in the string tension with coupling supports the conclusion that we are in a 
confining phase of the theory. Taking the Sommer scale to be ro = 0.5 fm, we find that the inverse 
lattice spacing ranges between 1.3 GeV at j8 = 2.3 and 2.1 GeV at j8 = 2.4. 
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Table 1: Static quark potential fit parameters for mf = 0.02 and = 16. 



j3 t range r range ro 

23 4^8 73^6 3.339(23) 

2.3533 5-9 ^3-6 4.379(80) 

2.4 5-10 V3-6 5.306(68) 



Vo « a 

0.501(18) 0.161(23) 0.134(4) 

0.569(23) 0.240(23) 0.074(4) 

0.539(11) 0.205(15) 0.051(2) 



□ P = 2.3533 




Figure 3: R as a function of time for mf = 0.02 and 
Ls = 16. 




2.3 2.32 2.34 2.36 2.38 2.4 



Figure 4: Residual mass as a function of coupling 
for mf = 0.02 and Lj = 16. 



3.2 Residual mass 

It is of crucial importance that we have an understanding of the residual mass (rrii-es) since 
its magnitude determines our proximity to the SUSY point in this lattice formulation. The 
dependence of the residual mass may be parameterized by the theoretically motivated formula : 

mres ~ ao— ha2^— , (3.3) 

where p(0) appearing in the second term (i.e. the dislocation term) represents the density of zero 
eigenvalues of the fifth dimension transfer matrix Hamiltonian. 

The residual mass was determined from a ratio R{t), given by the coupling of the pion ^ to the 
mid-point pseudo-scalar density divided by its coupling to the boundary (see [ ]T0[ ] for details). At 
large times this ratio of correlation functions is expected to tend asymptotically toward rrires', plots 
of this ratio for = 16 lattices are shown in Fig. (^. We determine the residual mass by fitting 
R{t) with a constant over the plateau region. Fig. (^ shows a plot of the extracted values of m^.^ 
as a function of the coupling. We find that the residual mass is roughly 5-10 times that of the input 
gluino mass. Furthermore, the strong coupling dependence of mres suggests that the dislocation 



term appearing in Eq. (3.3) dominates the residual mass. In order to reduce the residual mass, 



simulations at weaker couplings and larger Ls values are currently underway. 

3.3 Gluino condensate 

The gluino condensate was measured using a stochastic estimator with a single hit. We perform 



Here, we refer to the flavor non-singlet pseudo-scalar associated with a partially quenched theory. 
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Figure 5: Fits of the gluino condensate as a function of m/ (left; extrapolated values are indicated by x) 
and as a function of Ls (right; extrapolation error bars are indicated by dashed lines) for j3 = 2.3. 



chiral limit extrapolations of the gluino condensate using two different limit orders following 
First, we perform a linear mj ^0 extrapolation of the gluino condensate at fixed Lg, followed by 
an Li — > oo extrapolation of the mj = result using the best available, theoretically motivated fit 
formula: 



co + ci——, (3.4) 



which may be derived from the fifth dimension transfer matrix formalism. Note that the dislocation 
contribution to mres which appears in Eq. (3.3) is absent in Eq. (^^. This may be understood by 



observing that the chiral condensate is dominated by contributions from UV modes, whereas the 



dislocation term appearing in nires may be attributed to purely low energy phenomena |12p . 

Following the double extrapolation procedure outlined above we obtain an unrenormalized 
value of 0.003087(159) (x^/d.o.f. = 3.3) for the gluino condensate in the chiral limit at finite 
lattice spacing. Chiral extrapolations have been performed by reversing the order of limits and yield 
consistent results with comparable error bars. Plots of the m f and L., fits for these extrapolations are 
shown in Fig. (^. We have performed additional fits using other, phenomenologically motivated 
fit formulae to describe the dependence of the gluino condensate (e.g. Eq. ( |3.4| ), without the Lj^ 



prefactor in the second term). Such fits yield a 20% shift in the gluino condensate as compared to 



the value obtained with Eq. (3.4). 



3.4 Spectrum 

The low energy spectrum of ^ = 1 SYM is believed to consist of supermultiplets which may 
involve glue-glue, glue-gluino as well as gluino-gluino bound states. Although the states within 
a given multiplet are degenerate, at finite gluino mass (e.g. 7^ oo,mj 7^ 0) one expects mass 
splittings which are to leading order linear in the gluino mass. The mass splittings have been 
calculated using a variety of effective theories [13, |l^], however such calculations are unreliable 



since there is no separation of scales and therefore no small expansion parameter. 
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Table 2: Pseudo-scalar fit parameters for nif — 0.02 and Lj = 16. 



i8 


t range 


connected 
"'-ps 


t range 


Amps 


2.3 


9-23 


0.8701(3) 


2-4 


0.0180(37) 


2.3533 


9-23 


0.8144(6) 


2-4 


0.0176(52) 


2.4 


9-23 


0.7367(9) 


2-4 


0.0230(51) 



For the scalar and pseudo-scalar gluino-gluino and the glue-gluino composite states, we con- 
sider the interpolating fields 

il,(x,0 = TrA(x,Or,-A(x,0 , a{x,t) = TrF^y{x,t)L^yX{x,t) , (3.5) 

where F,- = {1,75} for the scalar (s) and pseudo-scalar (ps), X{x,t) is an appropriate interpolating 
field for the gluino, Jl^y = ^[7^, 7v] and F^y represents an interpolating field for the field strength 
tensor (e.g. a clover-leaf shaped product of link matrices). The lowest energy states created by the 
former operator correspond to the /o and rj' respectively in QCD, whereas there is no QCD ana- 
logue for the latter operator. Upon "integrating out" the fermion degrees of freedom, the correlation 
function for the scalar and pseudo-scalar operators involve a difference between two contributions: 
a "connected" part C,(f) and a "disconnected" part Di{t). As is the case with the rj' in QCD, the dis- 
connected contribution is numerically extremely difficult to evaluate exactly. We choose to instead 
use a stochastic estimator to approximate the correlator following the techniques of []T^]. 

The connected and disconnected correlators where measured using random Z2 volume and 
wall sources respectively (a single hit for the connected part and five hits for the disconnected 
part); to improve statistics we project onto the zero momentum state by averaging the result over 
all of space. In order to extract the mass of the pseudo-scalar, we study the ratio of disconnected 
and connected correlation functions 

Dpsit) 



Cpsit) 



2 — de 



(3.6) 



where Amps = nips ~ m'^^""'"^"''^ . Plots of this ratio at several different couplings are shown in Fig. 



(^). The linearity of these plots appear consistent with Ampst <^ 1, presumably due to the presence 
of a large residual mass. Assuming that this is the case, we may expand Eq. ( ^!^ ) to leading 
order in Ampst and then perform a linear fit to obtain a value for the mass difference Amps. Fig- (0) 
shows effective mass plots for the connected part of the pseudo-scalar correlation function at several 
different couplings from which ^^omecfed j^^y extracted. With Amps and /7jCon«ertet; determined, 
we may finally extract the pseudo-scalar mass mps. The results of these fits are provided in Table 
(^). While at present there are insufficient statistics to differentiate Amps between j8 runs, it is 
none-the-less evident that for each coupling mps is dominated by its contribution from m'^p""'^'^"^'^ . 
A complete analysis of the pseudo-scalar, scalar and their fermionic superpartner at smaller residual 
masses is currently underway and results will appear in a forthcoming publication. 
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Figure 6: Ratio of connected and disconnected Figure 7: Effective mass plot of the connected 
pseudo-scalar correlators as a function of time for pseudo-scalar correlator as a function of time for 
nif — 0.02 and Lj = 16. nif = 0.02 and L, = 16. 
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